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We study a one-dimensional lattice random walk with an absorbing boundary at the origin 
and a movable partial reflector. On encountering the reflector, at site x, the walker is reflected 
(with probability r) to x — 1 and the reflector is simultaneously pushed to x + 1. Iteration of the 
transition matrix, and asymptotic analysis of the probability generating function show that the 
' critical exponent 5 governing the survival probability varies continuously between 1/2 and 1 as r 

varies between and 1. Our study suggests a mechanism for nonuniversal kinetic critical behavior, 
i i . observed in models with an infinite number of absorbing configurations. 
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D 1 

Random walks find application in virtually every area of physics [1-5]. Random walks in the presence of traps or 
+Jj \ absorbing boundaries, and/or reflectors, are much studied as models of exciton recombination [6], diffusion- limited 

reactions [5] , and polymer-surface interactions [7] . 
c/3 \ In this Letter we study an unbiased random walk with an absorbing boundary at the origin, and a movable, partial 
reflector. At each encounter between walker and reflector, the latter moves one step to the right (i.e., to a site not 
yet visited by the walker), while the walker is reflected to its previous position with probability r. Thus the reflector 
hampers the advance of the walker into new territory, but does not otherwise influence its motion. (Note that in the 
limit r = the reflector has no effect, but it does mark the span of the walk.) We are primarily interested in the effect 
£^ ' of the reflector on the asymptotic scaling properties of the walk. 

The process defined above admits various physical interpretations. One is in terms of adlayer growth, with deposition 
and evaporation at the edge of the adlayer (but not in the bulk). The substrate adsorption sites are originally in 
a "non-activated" state, with a low sticking probability for incident molecules, but after a first adsorption attempt 
(at a step edge) the substrate site becomes activated, with a high sticking probability. Biological interpretations are 
also possible, e.g., of the advance of a bacterial colony in a growth medium, with a preliminary contact facilitating 
expansion into new regions, or, similarly, the spread of a political viewpoint in an initially skeptical population. Our 
results are of interest as an example of nonuniversality (a continuously-variable critical exponent) in a model that 
allows an exact asymptotic analysis. 

A further motivation for our study is provided by absorbing-state phase transitions, an area of great current interest 
in nonequilibrium statistical physics [8,9], in which a many-particle system, such as directed percolation, becomes 
trapped in a configuration allowing no further evolution. Continuous transitions to an absorbing state have been 
invoked in models of epidemics [10], catalysis [11], and the transition to turbulence [12,13], to cite but a few examples. 
The connection between absorbing-state phase transitions and random walks lies in compact directed percolation 
I ■ (CDP) [14]. CDP is defined on a discrete space-time (x,t) with time-slices corresponding to diagonals of a square 
lattice, so that the neighbors of site (x,t) at the previous time are (x—l,t—l) and (x+l,t— 1). Each site is either 
occupied or vacant. If (x, t — 1) has n occupied neighbors, then (x, t) is vacant (occupied) with probability 1 if n — 
(2), and is occupied with probability p if n = 1. CDP exhibits an absorbing-state phase transition at p = 1/2. 
Consider an initial state of a single occupied site in an otherwise empty lattice. The boundaries between the occupied 
region (descended from the initial seed particle) and the outer vacant regions follow simple random walks, which are 
unbiased if p = 1/2. Thus the length X(t) of the occupied region, being the distance between two random walks, is 
also a random walk, with an absorbing boundary at the origin. For p < 1/2 (> 1/2) X{t) is attracted to (repelled by) 
the origin. For p = 1/2, X(t) is unbiased, and well known results for random walks [1,4,15] imply that the survival 
probability decays for long times ~ t~ s , with 6 = 1/2. 

The motivation for introducing mobile reflectors in CDP (and thus in the simple random walk studied here) arises 
from the puzzling behavior of models that can become trapped in one of an infinite number of absorbing configurations 
(INAC) [16-18]. Anomalies in critical spreading for INAC, such as continuously variable critical exponents, have been 
traced to a long memory in the dynamics of the order parameter, p, due to coupling to an auxiliary field that remains 
frozen in regions where p — [18,19]. INAC appears to be particularly relevant to the transition to spatio-temporal 
chaos, as shown in a recent study of a coupled-map lattice with 'laminar' and 'turbulent' states, which revealed 



I 



> 

o 



o 
o 



c3 



l 



continuously- variable spreading exponents [20]. 

Grassberger, Chate and Rousseau (GCR) [21] proposed that spreading in systems with INAC could be understood 
by studying a model with a unique absorbing configuration, but in which the spread of activity to previously inactive 
regions is hampered (or facilitated). Our model represents a two-fold simplification of the GCR model: first, the 
appearance of inactive sites within a string of active ones is prohibited (analogous to going from DP to the more 
restrictive CDP process); second, we study a single random walker rather than the pair needed to describe CDP 
[23]. We believe, nonetheless, that our model captures the essential physics underlying anomalous critical behavior in 
models with INAC. 

We study an unbiased, discrete-time random walk on the nonnegative integers, x t = 0, 1, 2, with x = absorbing. 
Initially the walker is at xq = 1 and the reflector, whose position we denote by Rt, is at Rq~2. Each time the walker 
steps to the site occupied by the reflector, it is reflected back to R — 1 with probability r (and remains at R with 
probability 1 — r), while the reflector moves (with probability 1) to R+ 1. Evidently the process x t is non-Mar kovian, 
since the transition probability into a given site depends on whether it has been visited before. We can transform the 
model to a Markov process by enlarging the state space to include the reflector position; it is convenient to introduce 
the variable y t — Rt — l for this purpose. The process (x, y) is restricted to the wedge between x — (absorbing) and 
x = y, with transitions from y to y + 1 allowed only from the diagonal x — y. At this point it is useful to include a 
further generalization of our model, by assigning a probability p' for the walker to jump to the right (and q' = 1 —p' 
to jump to the left) when on the diagonal; these transition probabilities are summarized in Fig. 1. For p' < 1/2 
the walker experiences an additional impediment to visiting new territory, while for p'(l — r) > 1/2 the 'reflector' 
effectively becomes an accelerator, drawing the walker forward to a previously unvisited site. 

The nonzero transition probabilities for the Markov chain are: 

W[(x,y) - (x±l,y)] = l/2, x = l,...,y-l (1) 
W[(y,y) -> {y,y + l)]=p'r 
W[(y,y) ^ (y+l,y+l)}= p'(l-r) 

W[(y,y) - (y-l,y)] = l-p' . 

The probability P(x,y,t) follows the master equation 

P{x,y,t+l) = ±P(x-l,y,t) + ±P(x+l,y,t), x = 1, 2, . . . , y-2; y > 3 , (2) 

with P(Q,y,t) = 0, representing the absorbing boundary at x = 0. Letting D(y,t) = P(y,y,t), the boundary 
conditions along the diagonal are: 

P(y-l,y,t+l)= ±P(y-2,y,t)+p'rD(y-l,t) + (l-p')D(y,t) , (3) 
D(y,t+1) = \P{y-l,y,t)+p'{l-r)D{y-l,t) , 

for y > 3. For y = 1, 2 the evolution equations depend on the initial condition; here P(x,y,0) = 5 x ,i5 y ,i [24]. 
Define the generating function P(x, y, z) — Y^tLv z * P{ x i U> *) ( an d similarly for D 7 etc.). P satisfies 

z- 1 P(x,y) = ±P(x-l,y) + ±P(x+l,y), x = 1, 2, . . . , y - 2; y > 3 . (4) 
(we drop the argument z for brevity), subject to the boundary conditions 

P(0,y) = 0, (5) 
z^Piy-l.y) = ±P(y-2,y)+p'rD(y-l) + (l-p')D(y) , (6) 
z-'Diy) = lP(y-l,y)+p'(l-r)D(y-l) . (7) 

The solution of Eqs. (4) and (5) is 



P(x,y)=C(y)(X x + -X x _); \ ± = z' 1 ± ^ z^ - I , (8) 

with C{y) yet to be determined. Noting that A + = XZ 1 = X, we define X(x) = X x — X~ x ; thus P(x,y) = C(y)X(x). 
Note also the recurrence relation for integer x: X(x— 1) — 2z~ 1 X(x) + A(x+1) = 0. 
^From Eq. (7) we find 

v 

b{y)=w y - 2 b{2) + \zu y Y J UJ- y ' P{y'-l,y') , uj=p'(l-r)z. (9) 

y'=3 
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It follows that D(y) — uD(y—l) = ^C(y)X(y — 1). Therefore, subtracting from Eq. (6) uj times the corresponding 
equation for y— 1, and using the recursion relation for X(x), we find 

C(y) , rA(y-2) + (l-r)A(y-l) 

~P 73 1TT7 TT—TT?— TT > ( 10 ) 



C(y-1) y {p'-\)X{y-l) + \X{y+l) 

which, noting that C(2) = P(l, 2)/A(l), provides the complete solution for the probability generating function. 

To determine the survival probability S(t) = J2"^Lo J2x=o P{ x i f> *) as * ~~ * °°' we ana ly ze the singular behavior of 
S'(z) as z — > 1. In this limit A = 1 + \/2e + 0(e), where e = 1 — z, and we find the dominant term to be [22] 

oo 

]T(7( y )[A(y/2)] 2 , (11) 



\/2(l-z) 



which converges for z < 1. Note that A(n) ~ 2n\/2e for ne <C 1, while in the opposite limit A(n) ~ A". For m< 1 
we may thus use the z — ► 1 limiting expression for C(n), that is, 

^(1,2) r(2 + W) r( y -r) 

C ° (2/) ~2 v /2(l-z) l'(2-r) 1-( W +WJ - (12) 
We then approximate the summand in Eq. (11) by 



2n 2 C* (n)e, ne < 1 

27V 2 C*o(iV)ea n - JV A rl ~ Ar , ne > 1 



1 "*T2f<„(AT\^s.n-N\n-N - 1 ; " >; 



where N — 1/y/e rounded to an integer, and 

, r + (l-r)A 

a_P (p'-l/2)A+(l/2)A3 ' lMj 

Since S and F(z) = X)^Lo /« have the same radius of convergence and linin^oo (7(n)[A(n/2)] //„ is finite, an 
Abelian theorem implies that the singular behavior of S is the same as that of F(z)/y/e [4]. Write F(z)/\fe = 
2c 1 / 2 J2n=o n2 Co(n) + 2e 1 / 2 N 2 Co{N) Y^=N+i(' Ka ) n ~ N ■ We approximate the first sum by an integral; using (%(n) ~ 
e -i/2 n -r-i/p we ^ nen g n( j ^ ne g rs ^ term <~ e( r + 1 /p ^ 3 )/ 2 . The same holds for the second term when we note 
that 1 - aX ~ e 1 / 2 and Co{N) - e -5( 1+r+1 /f'). Thus 5 ~ (1 - z) 5 " 1 (plus less singular terms) as z -> 1, where 
<5 = (r + 1/p' — l)/2. Finally, noting that for large n, the coefficient of z™ in (1 — z) 8 ^ 1 is n~ s /T{1 — S), we have that 
the survival probability decays asymptotically ~ t~ s ; in particular, for p' = 1/2 we have ~ t _ ( 1+r '/ 2 . Note also 
that the larger p' is, the slower the decay of the survival probability, as expected. 

While the amplitude of the leading term in S(t) can be evaluated exactly for p' = 1/2 and r = or 1 [22], we 
have not found a closed-form expression in the general case. Numerical evaluation of the exact expression for S does, 
however, lead to an amplitude in perfect agreement with that found in the transition matrix analysis. 

We turn now to results obtained via numerical iteration of the transition matrix defined in Eq. (2). We studied 
the process to a maximum time of up to 3 x 10 5 , quite sufficient to observe asymptotic behavior. In Fig. 2 we show 
the survival probability for p' = 1/2 and various values of r. In each case there is a power- law decay, with the decay 
exponent S varying continuously with r. (The plateaus in S(t) at short times reflect that for r = the walker can only 
reach the origin on odd-numbered steps. Similarly, for r = 1 the walker cannot reach the origin on the third step.) 
We verify that the mean position (xt) ~ t 1 / 2 (and that (x 2 ) ~ t; the averages are over surviving walks) independent 
of r. The amplitudes for the displacement and its second moment do decrease with increasing r. 

To determine the decay exponent precisely, we study the local slope S(t), given by a least-squares linear fit to the 
In 5 data for a set of 11 equally-spaced values (increments of 0.05) of Int. Plots of S(t) vs. t~ x lead to asymptotic 
(t — > 00) values confirming 6 = (l + r)/2 to better than one part in 2000. It is possible to improve the analysis further 
by determining the leading correction to scaling in an expansion of the form: S(t) ~ At~ 5 [l — Bt~ Al + •••]. For r = 

we find A = 0.7979, B = 0.78 and A 1 = 1.00, in accord with the exact expansion P(t) ~ y/2/(-Kt)[l - (3/4)^ H ]. 

For r = 1 our data yield A = 1.0000 (in agreement with the exact result), B = 1.332 and Ai = 1.00. In these two 
cases, plotting 8(t) versus t~ Al = t^ 1 we obtain 6 = 0.5000. and 0.999985, respectively. For r in the interval [0.1,0.8] 
we find a much smaller correction to scaling exponent, Ai ~ 0.54, while for r = 0.9, Ai ~ 0.62. Plotting S(t) versus 
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i~ Al (see Fig. 3), we obtain exponents that agree with the theoretical value to better than one part in 2 x 10 4 . The 
attractive simplification, Ai = 1/2 for r / 0, 1, yields essentially the same values for S, and is in fact to be expected, 
given that generic corrections to S should be cx (1 — z ) <5 - 1 +2 m (pi = 1,2,...). 

The case of compact directed percolation is somewhat more complicated than the random walk problem analyzed 
above, as we must now keep track of three variables, viz., the distance X between the two walkers, and the respective 
distances between the walkers and the associated reflectors [23]. Simulation results for CDP with movable partial 
reflectors again show 5 varying continuously with r, but over a somewhat wider range: for r = 1, 5 = 1.158(3). 

In one-dimensional spreading models (contact process [10], pair contact process [16], etc.) the survival and spread 
of activity may again be described in terms of the position of two points that bound the active region, as in CDP. But 
here the boundaries of the active region do not follow simple random walks: their dynamics involves a variable step 
size (due to gaps in the distribution of active sites), and a significant memory (the step-size distribution depends on 
the history of the previous step directions). It seems reasonable, nevertheless, to expect that varying the probability 
for advancing into virgin territory will change the scaling of the survival probability, just as observed here. A detailed 
investigation of this issue is a high priority for future work. 

In summary, we have uncovered the rather remarkable property of a continuously varying critical exponent governing 
the survival probability of a random walk with an absorbing boundary, in the presence of a movable partial reflector. 
Our results suggest a new approach to understanding nonuniversal spreading in models with an infinite number of 
absorbing configurations. 
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Figure Captions 

FIG. 1. Transition probabilities in the (x,y) plane. 

FIG. 2. Survival probability versus time for random walk with a movable partial reflector, p' — 1/2. Reflection 
probability (top to bottom) r = 0, 0.1, . . . , 1. 

FIG. 3. Local slope S(t) for p' = r = 1/2. The inset shows the same data plotted versus 1/t. 



5 




X 




2 4 6 8 10 12 

Int 



■0.740 



-0.742 



■0.744 



-0.746 



-0.748 



-0.750 




0.00 0.01 0.02 0.03 0.04 0.05 

f-0.53 



